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Abstract. Present day estimates of the Hubble constant based on Cepheids and on the cosmic 
microwave background radiation are uncertain by roughly 10% (on the conservative assumption 
that the universe may not be perfectly flat). Gravitational lens time delay measurements can 
produce estimates that are less uncertain, but only if a variety of major difficulties are overcome. 
These include a paucity of constraints on the lensing potential, the degeneracies associated with 
mass sheets and the central concentration of the lensing galaxy, multiple lenses, microlensing by 
stars, and the small variability amplitude typical of most quasars. To date only one lens meets all 
of these challenges. Several suffer only from the central concentration degeneracy, which may be 
lifted if one is willing to assume that systems with time delays are either like better constrained 
systems with non-variable sources, or alternatively, like nearby galaxies. 



1. Context 

Were we to choose at random a hundred members of the lAU and ask each of them to 
tell us the value of the Hubble constant and how it is measured, few if any would refer to 
gravitational lens time delay measurements. Most would instead refer to observations of 
Cepheid variables with HST or to observations of the CMB power spectrum with WMAP. 
These set the de facto standard against which time delay estimates must be evaluated. 

The Cepheids give a Hubble constant of 72 km/s/Mpc with a 10% uncertainty (Freed- 
man et al. 2001), a result that most of our randomly chosen astronomers would find 
relatively straightforward. But only a handful of them would be able to tell us how the 
CMB power spectrum yields a measurement of the Hubble constant. 

In his opening talk, David Spergel told us that three numbers are determined to high 
accuracy by the CMB power spectrum: the height of the first peak, the contrast between 
the heights of the even and odd peaks, and the distance between peaks. Ask anyone who 
calls himself a cosmologist how many free parameters his world model has and the number 
A'^ will be larger than three. Spergel's three numbers constrain combinations of those A^ 
parameters, but do not, in particular, tightly constrain the Hubble constant. One must 
supplement the CMB either with observations, or alternatively, with non-observational 
constraints. 

The effects of supplementing the observed CMB power spectrum with other observa- 
tions, and of adopting (or declining to adopt) non-observational constraints, are shown in 



table 1. All of the numbers shown are taken from the analysis by Tegmark et al. (2004b) 
They use the first year data from the Wilkinson Microwave Anisotropy Probe (Bennett 
et al. 2003), the second data release of the Sloan Digital Sky Survey (Abazajian et al. 
2004), the Tonry et al. (2003), data for high redshift Type la supernovae, and data from 
6 other CMB experiments: Boomerang, DASI, MAXIMA, VSA, CBI and ACBAR. 

f Present address: Center for Space Research 37-664G, 77 Massachusetts Avenue, Cambridge, 
MA 02139, USA 
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Table 1. The CMB power spectrum and the Hubble constant 
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We see from the first line of tabled that the value of the Hubble constant derived from 
the WMAPl data alone is a factor of 1.5 or 2.5 more uncertain than the Cepheid result, 
depending upon whether the universe is taken to be perfectly flat or only very-nearly flat. 
Combining the galaxy power spectrum as measured with SDSS2 (Tegmark et al. 2004a) 
reduces the uncertainties by factors of 3 and 2, respectively. In the very-nearly flat case, 
adding type la supernovae does little to change the uncertainty but shifts the value of Hq 
closer to the Cepheid value. In the perfectly flat case, adding other CMB measurements 
reduces the uncertainty in Hq by another factor of 1.3. Small deviations from flatness 
produce substantial changes in the Hubble constant, with the product hfl^^^ remaining 
roughly constant (Tegmark et al. 2004b). 

Suppose we were to show tabled to our randomly chosen lAU members and again ask 
them the value of the Hubble constant (and its uncertainty) . Some would argue that the 
results for the very-nearly flat case are so very-nearly flat that it is reasonable to adopt 
perfect flatness as a working model. We observe that | log^ltot \ < 0.10, when it might have 
been of order 100 (or perhaps 3 if one admits anthropic reasoning). But others would be 
reluctant to take perfect flatness for granted. The very sensitivity of the Hubble constant 
to that assumption would be an argument against adopting it. The present author is, 
himself, ambivalent. 



2. Ho from time delays circa 2003 



Refsdal (1964) first proposed measuring the Hubble constant using gravitationally 



lensed supernovae 15 years before the discovery of the first gravitationally lensed quasar. 
In the past two years the status of Hubble constant measurements from time delays has 
been reviewed by Courbin (2003) and by Kochanek & Schechter (2004) Figure taken 
from Courbin's review, summarizes his report. It should be noted that Courbin deliber- 
ately incorporated results representing diverse approaches to computing Hq. The average 
is one standard deviation smaller than the Cepheid result with a comparable uncertainty. 
It is also in agreement with the CMB results of tabled 

But if we were to make such a figure today, we might find it less reassuring. The Hq 
value for PG1115-I-080 in figureHis based on a model for the lensing galaxy that is far 
from the fundamental plane for ellipticals. The value for RXJ091 1+0554 is based on a 
faulty model (Schechter 2000) that the present author repudiates in M.4I And the values 
for HE1104-1805 are based on a time delay that |Ofek fc Maoz (2003) have shown to be 
too large by a factor of two, driving the points off the top of the plot. The uncertainty in 
the Hubble constant, as derived from the scatter among systems, would be considerably 
larger today than it is in figure Q 

In what follows we discuss a number of major difficulties associated with using time 
delay measurements to measure Hq. All but one of the systems in figure ^ are subject 
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Figure 1. Time delay measurements of tiie Hubble constant, circa 2003 (Courbin 2003). 



to one or more of these difficulties. Before delving into the details, we review the physics 
upon which the estimate of the Hubble constant is based. f 



3. Lensing fundamentals 

3.1. The 3 D's 

There is a sense in which time delay is the most fundamental manifestation of gravi- 
tational lensing. In the weak field limit a gravitational potential produces an effective 
index of refraction that increases the travel time for a photon. Fermat's principle requires 
that the photon travel along a path that is a minimum, a maximum or a saddlepoint 

f An exc ellent, although slightly outdate d pedagogical treatment of gravitational lensing can 
be found in |Narayan fc Bartelmann (1999)] 
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of the travel time. The photon must detour around the path it otherwise would have 
taken. A balance is struck between increased travel time associated with the detour (the 
geometric delay) and decreased travel time due to the shallower gravitational potential 
(the Shapiro delay), 



/ time \ / geometric\ / gravitational\ 
\delay/ \ delay / \ delay / 



(3.1) 



If the size of the lens is small compared to the path length, we can use the thin lens 
approximation. The time delay can be expressed in terms of a two-dimensional effective 
potential obtained by integrating the gravitational potential $3/5 along the line of sight, 

J observer ^ ^ L 

It depends upon the angular diameter distances to the lens and to the source, and 
Ds, and upon the distance from the lens to the source, D^s- With this simplification all 
of gravitational lensing boils down to taking derivatives of the time delay function. 



1 + zl dLds 

T = 



l{e-pr-^2Die) 



(3.3) 



where 9 measures the position on the sky at which a ray crosses the plane of the lens and 
/3 is the position on the sky of the source. Here we use the dimensionless variants of the 
angular diameter distances, dz,, ds, and dLS- The quantity — P is the amount by which a 
ray is deflected. The geometric part of the time delay varies as the square of deflection in 
the limit of small deflections, a straightforward consequence of the Pythagorean theorem. 

Differentiating r with respect to 6 and setting the gradient of r to zero allows us to 
solve for the minima, maxima and saddlepoints of the time delay. It gives what is often 
call the lens equation, 

e- P-V^/J2D = . (3.4) 

The deflection, 6 — [3, is equal to the gradient of the two dimensional potential. 

The magnification of an image is the ratio of the observed angular size to the true 
angular size of an object. We want something like the ratio of A9 to A/3, but both of 
these are vectors. Passing to derivatives, the quantity dO/dfi gives a magnification matrix 
that maps a vector at the position of the source, A/3, into a vector in the image plane, 
A0. This matrix cannot be calculated from equation 13.41 because the two dimensional 
potential is a function of position in the lens plane, not the source plane. But its inverse 
can be calculated by taking the derivative of equation 13.41 with respect to 9, giving the 
inverse magnification matrix, 

(3.5) 





This is readily inverted to give the magnification matrix, 89/8(3. Since the magnification 
matrix is symmetric, there is a choice of coordinates for which it is diagonal. Those 
diagonal elements are in general unequal, stretching or compressing the image more in 
one direction than the other. The magnification matrix might equally well be called the 
distortion matrix. 

It is often the case that the source and its images are very much smaller than the 
resolution of our observations. Gravitational lensing conserves surface brightness, but as 
the solid angle subtended by the image is usually different from that subtended by the 
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source, the observed flux is also different. The scalar magnification factor, /i, is given by 

dd 



^,= — . (3.6) 

op 

The observable consequences of lensing are mnemonized with three "D's," each cor- 
responding to a different derivative of the time delay function. The first derivative of 
the delay function gives a deflection. The second derivative gives a distortion. And the 
"zeroth" derivative is just the delay function itself. 

3.2. Modeling the lens 

The previous subsection illustrates the power of mathematics: an enormous amount of 
interesting astrophysics can be swept under the rug of a single function, the two di- 
mensional gravitational potential, ip^O). If we knew a priori we would long since have 
measured the Hubble constant. Instead we must construct models for our potentials, ad- 



justing them to fit whatever observations can be brought to bear. Williams & Saha (2000) 
have developed a method for modeling lenses that avoids making detailed assumptions 
about the shape of the lensing potential (see also SAHA & Williams' contribution to the 
present proceedings), at the price of non- uniqueness. The present author takes the view 
that we know something about galaxy potentials and therefore would do well to start by 
assuming lensing galaxies are much like nearby galaxies. 

A model which is both very simple and very useful is the singular isothermal sphere 
in the presence of an external tide, 

MO) = br + ^r^ cos2{(j) - (by) . (3.7) 

The leading (monopole) term is the projection of the isothermal's potential. In our case 
the strength of the isothermal is parameterized by 6, the radius of its Einstein ring. The 
tidal term is a quadrupole, and is parameterized by a dimensionless strength, 7, and an 
orientation (p^. Note that here we represent angular position on the sky, 0, in terms of 
polar coordinates r and cj), with r measured in radians. 

The singular isothermal sphere gives the flat rotation curves characteristic of nearby 
galaxies. A substantial fraction of quadruple gravitational lenses are the result of strong 
tides from a group or cluster of galaxies of which the lens is a member. In the absence of 
these tides the lenses would produce only double images. The intrinsic flattening of the 
lensing galaxy also produces a quadrupole term (which has a different radial dependence), 
but the tidal contribution is usually considerably larger. 

Our simple model does great injustice to the complexity of individual lenses, but 
complexity is a luxury we cannot afford. A parameterized model can have at most as 
many parameters as the available observable quantities and we are quite limited. The 
numbers of observables are summarized in table [21 Delay is a scalar, deflection is a two 
dimensional vector and the symmetric distortion matrix has three independent elements. 
We have subtracted one deflection from our available constraints because we do not know 
the position of the source. Likewise we subtract one distortion because we do not know 



Table 2. Constraints 



Delays: 1 x [^images — 2] 

Deflections: 2 x [^images — 1] 
Distortions: 3 x [^images — 1] 

(or for unresolved sources) 
magnifications: 1 x [^images — 1] 
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the size and shape of our source (taking it to be elhptical). By the same argument we 
subtract one delay because we do not know the absolute time at which some event used 
to calculate delays actually occurred. In the case of delays we subtract a second delay 
because we need it to calculate the Hubble constant. It cannot be used both to constrain 
our model and to give Hq. 

The situation is even worse if the images are not resolved by our telescopes. In such 
cases we measure only the fluxes of the images and get only one number from each 
distortion rather than three. Consider the case of a lensed quasar with two unresolved 
images. After discounting for the unknown properties of the source and reserving one 
delay for measuring Hq we have only three constraints for our model. Using even our 
very simple model we are guaranteed a perfect fit, with no additional degrees of freedom 
left to check goodness of fit. 

With no sanity check on our model, such doubles can be used only with extreme 
caution. Something might be amiss and we would have no way of knowing it. Moreover 
we cannot add additional parameters (allowing, for example, for the shape of the galaxy) 
because we have exhausted our supply of constraints. Four of the nine systems in figure^ 
are just such doubles: B1600+43, HE2149-2745, SBS1520+53 and HE1104-1805. 

The situation is slightly improved for four image systems. It is likewise slightly im- 
proved if there are multiple sources, as is sometimes the case for radio AGN. In con- 
structing a list of major difficulties faced in measuring Hq, the first would be the paucity 
of constraints. 



4. Problem lenses 

4.1. Kochanek's time delay formulation 

The system PG1115+080 illustrates another of the major difficulties in time delay mea- 
surements of Hq. The mathematics is simpler if we diagonalize the magnification matrix 
and introduce two dimensionless quantities, the convergence n and the shear 7. 








(4.1) 



l-K+7 / 

The convergence is the two-dimensional Laplacian of the two-dimensional potential, 

K-iv^V , (4.2) 

and is therefore proportional to the surface mass density. 

For zero shear, a convergence of unity gives infinite magnification in both directions. 
For zero convergence (which means zero surface density), a shear of unity gives infinite 
stretching in one direction. For a point mass there is unit shear everywhere on its Einstein 
ring. 

KOCHANEK (2002; present proceedings) has shown that with these definitions we can 
write an expression for the difference in time delay between two images in a circularly 
symmetric system that depends to first order only upon the positions of the images with 
respect to the lens and the convergence averaged over the annulus bracketted by the two 
images. 



1 
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Figure 2. The quadruple lens PG1115+080, observed with the Baade 6.5-m telescope. The 
small group of galaxies to the upper right lies very close to the position predicted from the 
positions of the four quasar images. It lies at the same redshift as the lensing galaxy. 

For power law density profiles, p ^ r"'', we have 

(-)-^ . (4.4) 

Notice that ii rj — 2, as would hold for an isothermal sphere, the average value of the 
convergence is one half. As 77 ^ 3, approaching the central concentration for a point 
mass, (k) — + and the predicted time delay is twice as long as in the isothermal case. 

4.2. The mass sheet degeneracy 

Kochanek's parameterization illuminates two different difficulties in measuring Hq from 
time delays. Figure |21 shows PG1115+080 and its immediate surroundings. A fit of our 
simple model to the image positions indicates that the source of the tide lies along a line 
passing through the lensing galaxy and the upper right corner of the figure. At the cost 
of one additional parameter, we can substitute a second isothermal sphere for the tidal 
shear. This very much improves the fit, and the derived position is coincident with the 
small group of galaxies to the upper right of the lensed system. The inferred velocity 
dispersion for the second isothermal sphere is 375 km/s (Schechter et al. 1997), typical 
of a small group. The observation of a group at the position "predicted" by the simple 
model inspires confidence in it. 

But if the group is isothermal, it has a convergence associated with it - in this case ap- 
proximately 10%. This decreases the predicted time delay by 10%. Small numbers make 
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Table 3. Exponents for power law density profiles: p ^ r 



lensing galaxy 77 reference 


JVAS0218+357 IMt^^^ 
Q0957+561 1.84 
MG1131+0456 2.40to;2 
PMN1632-0033 1.9lloo2 
MG1654+1346 1.90tgjf 
CLASS1933+503 1.86l^iJ 


Wucknitz, Biggs, & Browne (2004)| 


Barkana ct al. (1999) 


Chen, Kochanek, & Hewitt (1995)| 


Winn, Rusin, & Kochanek (2004)| 


Kochanek (1995)| 


Cohn, Kochanek, McLeod, & Keeton (2001)| 


< > 1.98 



measurement of the actual velocity dispersion difficult, but it is consistent (Tonry 1998) 
with the tidal estimate. Many of the quadruple lenses live in groups and clusters with 
convergences at the position of the lens of 10, 20 and even 30%. The corrections are sub- 
stantial, and we must wonder whether we have got the cluster convergence correct. Two 
of the systems in figure ^ Q0957+561 and RXJ091 1+0554, have inferred convergences 
in the 20 - 30% range. This is a manifestation of a difficulty more generally known as the 
mass sheet degeneracy (Gorenstein, Shapiro, & Falco 1988; Saha 2000). There is no way 
of knowing (from deflections and distortions) whether or not a mass sheet is present, but 
it affects the predicted time delays by changing the mean convergence. 



4.3. The central concentration degeneracy 

Eauations l4.3l and l4.4l tell us that the interpretation of the time delay for PG1115+080 de- 
pends crucially upon the internal structure of the galaxy. An error in the radial exponent 
of plus or minus 0.1 produces a 10% change in the predicted time delay as compared 
to an isothermal model. Fortunately we can test the isothermality hypothesis. For an 
isothermal sphere, the angular radius of the Einstein ring, 6, is a directly proportional to 
the square of the velocity dispersion, a^: 

b^^in^ (4.5) 
as 

With an Einstein ring radius b = 1'.'15, a lens redshift = 0.31 and a source redshift 
zs = 1.71 we predict cs/s = 232 km/s. This must be reduced by a factor -^(1 — Hgroup) 
which amounts to a 5% correction in the present case. We can test our prediction, but 
the measurement is a difficult one. The lensing galaxy is crowded by four very much 



brighter images. Tonry (1998) measures a velocity dispersion of 281 ±25 km/s, very much 



larger than predicted. Treu & Koopmans (2002) use Tonry's measurement in modeling 



PG1115-)-080 and derive a value for Hq very much larger than under the isothermal 
hypothesis. Their value is the one plotted in figure ^ 

Why are we so reluctant to abandon the isothermal hypothesis for PG1115+080? 
First, because velocity dispersion estimates from equation 14. 51 for an ensemble of lensing 
galaxies are consistent with the fundamental plane relation for non- lensing ellipticals. 
PG1115+080 is in no way unusual (Kochanek et al. 2000), but would be if we adopted 
the direct measurement. 

A second argument is that the lenses for which we can measure the radial exponent are 
very nearly isothermal. Lenses with multiple sources, rings or central images all break 
the central concentration degeneracy and permit measurement of the radial exponent rj. 
The results for six systems are shown in table 13 (see also Wayth's contribution to the 
present proceedings). While there are differences in the way rj was calculated for each 
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Figure 3. The quadruple system HE0230-2130, observed with the Baade 6.5-m telescope 
through g' (left) and i' (right) filters. The two red objects are lensing galaxies. The image 
between the two galaxies is a saddlepoint of the time delay function. 



of these systems, the results are consistent with isothermal. Unfortunately only one of 
these systems, JVAS0218+357, is also a system that has a measured time delay. 

A third argument is that a great many nearby galaxies have been studied, and their 
potentials are consistent with isothermals. This is nicely shown in a figure published by 
Romanowsky & Kochanek (1999) For a sample of twenty bright ellipticals, the circular 



velocity inferred from the velocity dispersion declines only slightly over a factor of 30 in 
radius. The corresponding decline in PG1115+080, computed from its measured central 
dispersion and the Einstein ring radius, is very much larger. It is nonetheless true that 
the potentials for nearby ellipticals are only very-nearly isothermal and not perfectly 
isothermal. The central concentration degeneracy qualifies as another major difficulty 
associated with time delay estimates of Hq. 



4.4. Multiple lenses 

Another problem associated with modeling lensing galaxies is illustrated by the system 
HE0230-2130. There is not one lensing galaxy, but two, as seen in figure |31 The second, 
fainter galaxy provides the quadrupole moment that makes this a 4- image system and not 
a 2-image system. When we see two galaxies separated by only a few kpc, we must wonder 
how the non-baryonic matter is distributed. Are the halos for the two galaxies distinct? 
Or have they merged? How many components should we model? And where should 



their centers be? In such a case something with elements of the Williams & Saha (2000) 



form-free approach might be preferable to a straightforward parameterized model. 



The system with the best time delays is CL ASS 1 608+656. Fassnacht, Xanthopoulos, Koopmans, & Rusin (2002) 



have measured four beautiful radio lightcurves. All four time series faithfully reproduce 
the many bumps and wiggles. The multiple delays in this case are good enough to provide 
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Figure 4. The quadruple system CLASS1608+656, observed with the Hubble Space Telescope 
through the FSlAW filter. The lensing galaxies are marked Gl and G2. A dust lane can be seen 
encircling them. 



meaningful additional constraints to the deflections and distortions. But modeling the 
system is not straightforward. As with HE0230-2130, there are two lensing galaxies (fig- 
ure^J) that appear to be interacting. A dust lane encircles them both. Again the question 
of how the dark matter might be distributed looms as crucial. Koopmans et al. (2003) 
have worked exceedingly hard to constrain this system, measuring positions for the im- 
ages, the shape of the ring, and measuring the velocity dispersion of the more massive 
lens. We can nonetheless imagine a devil's advocate coming up with a plausible dark 
matter distribution uncorrelated with the observed galaxies that gives a very different 
value for the Hubble constant. 

The problem of multiple lenses can be serious even when one of the galaxies is very 
much smaller than the principal lens. Consider the case of RXJ091 1+0554, shown in 
figureEl where the primary lensing galaxy has a faint companion. Given its faintness and 
our aversion to adding additional parameters, we might be tempted to ignore it, as did 



Schechter (2000) in the model used by Courbin to construct figure ^ But allowing for a 
mass at the position of this dwarf companion changes the predicted time delay by 10%. 
Though the smaller galaxy is a factor of 10 fainter than the primary lensing galaxy, its 
effect is to move the center of mass closer to the midpoint of the images, decreasing the 
differences in path length. 

There is an irony here in that the multiple lenses work to our advantage in producing 
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Figure 5. The quadruple system RXJ0911+0554, observed by the CASTLES consortium with 
the Hubble Space Telescope through the F160W filter. The primary lensing galaxy Gl has a 
dwarf satellite G2. Allowing for the mass of the dwarf changes the derived value of Ho by 10%. 

(or adding to) the quadrupole moments that give us 4 images rather than 2. But at the 
same time they make modeling considerably more difficult. 

4.5. Microlensing 

In his contribution to the present proceedings, Lutz WiSOTZKl presents spectroscopic 
evidence for microlensing of quasars. This has important consequences (unfortunately 
negative) for time delay estimates of Hq . It means that there are uncorrelated variations 
in the fluxes of images, especially in the optical continuum. The timescale for these 
fluctuations depends upon the size of the Einstein ring of the microlenses and upon the 
relative' velocity of the quasar and the lensing galaxy. Taking the lenses to be a solar mass 
and the relevant velocities to be of order 300 km/s gives timescales of order 10 - 30 years. 
This ought not to matter for quasars lensed by galaxies, whose time delays are at most 
one year. But there have been repeated instances of uncorrelated variations in quasar 
lightcurves on considerably shorter timescales (e.g. Burud et al. 2000), in some cases 
causing considerable dispute about which points ought and ought not to be included in 
a time delay measurement and about the nature of the fluctuations. 

Seven years ago the present author prevailed upon Andrzej Udalski to monitor HE1104- 
1805 for a time delay with the OGLE telescope (Udalski, Kubiak, & Szymanski 1997). 
The OGLE group was at that time focussing exclusively on LMC and galactic bulge 



12 Paul L. Schcchtcr 

600 800 1000 1200 1400 1600 1800 2000 2200 

n \ \ \ \ \ \ \ \ 




600 800 1000 1200 1400 1600 1800 2000 2200 

HJD^ - 2450000.0 



Figure 6. Lightcurves for the two components of HE1104-1805, taken with the OGLE 1.3-m 
telescope (Wyrzykowski et al. 2003). Note that while the measurement uncertainties are larger 
for the fainter, B image, the scatter is larger for the brighter, A image. 



microlensing, but HE1104-1805 is at a very different right ascension. Three years of superb 
OGLE data were obtained in which no believable correlated variations were observed. But 
there were appreciable uncorrelated variations, most of which were attributable to the 
brighter of the two images. After three years we admitted defeat and wrote a paper about 
microlensing (Schechter et al. 2003) which was, after all, OGLE's primary mission. The 
OGLE telescope was out of commission for upgrades during much of the following year; a 
few data points were obtained but were never reduced. But by then Ofek & Maoz (2003) 



had begun monitoring HEl 104-1805 at the Wise Observatory. They used our data in 
conjunction with theirs to give a time delay that differed by a factor of two from the one 



used in figure ^ After Ofek and Maoz published their delay, Wyrzykowski et al. (2003) 
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reduced the remaining OGLE data, confirniing their value and giving the hghtcurves in 
figure 

Notice that the error bars are smaller for the A image, yet the lightcurve is smoother 
for the B image. Image A appears to vary with an amplitude of 0.06 magnitudes on a 
timescale of one week. The interpretation of these fluctuations is somewhat speculative, 
but there is no question that microlensing represents an another difficulty for time delay 
measurements of Hq . 

4.6. The myth of quasar variability 

One last difficulty associated with estimating Hq is implicit in the fact that only 10% 
of the hundred or so known lensed quasars appear in figure When we think of quasar 
variability we think of 3C273 or 3C279 and variations of a magnitude or more. But 
those quasars are well known precisely because of their variability. Most quasars are 
considerably more boring. Systematic studies of quasar variability at optical wavelengths 
consistently give rms fiuctuations of order 10 - 15% on proper timescales of one year 
(Cristiani et al. 1996; Vanden Berk et al. 2004). 

The good news is that a few quasars vary more than this. The bad news is that most 
vary even less than this. The statistics may be somewhat better for radio quasars, but then 
radio loud systems constitute only 10% of all quasars. While the small amplitude of quasar 
variations may give graduate students (mostly those working at radio wavelengths) an 
opportunity to show their skill in beating down their observational errors, flat lightcurves 
don't lead to offers of prestigious postdocs. It should be noted that there is considerable 
room for improvement in the optical lightcurves, for which the observational errors are 
much larger than photon statistics would imply. Our list of major difficulties in measuring 
Hq has grown to include the following: 

• the paucity of constraints, 

• the mass sheet degeneracy, 

• the central concentration degeneracy, 

• multiple lenses, 

• microlensing and 

• small variability amplitudes. 



5. A golden lens 

The difficulties we have described are sufficiently common that almost every lensed 
quasar is subject to at least one of them. In general (though there are important excep- 
tions) the double systems suffer from a paucity of constraints. Many quadruples suffer 
either from the mass sheet degeneracy (due to a nearby cluster) or from multiple lensing 
galaxies. Only a few systems (mostly radio loud quasars) have multiple sources permit- 
ting unambiguous determination of a radial density exponent. Are there any lenses which 
are not obviously subject to one or another of these difficulties? 

A decade ago we imagined finding a "perfect" (Press 1996) or "golden" lens (Williams 
& Schechter 1997) that would permit determination of Hq with small uncertainty. A 
lens that was not subject to any of the above difficulties might reasonably qualify for 
certification as a golden lens. 

The system JVAS0218-I-357 (Patnaik et al. 1992) is such a system. It is a radio double 
that, at VLBI resolution, has a core and a blob (Biggs et al. 2003) of the sort that radio 
astronomers fancifully call jets. The double source breaks the concentration degeneracy, 
and implies an isothermal potential. _Wucknitz, Biggs, fc Browne (2004) have modeled 



the radio data, finding the radial exponent and the center of symmetry of the lensing 
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Figure 7. The shading of the pixels indicates relative likelihood for the optical position of 
the galaxy lensing JVAS0218+357 (York et al. 2004). The solid line is drawn from the B im- 
age toward the A image, which lies beyond the plot boundary. The diagonal dashes indicate 
loci of constant Hubble constant, starting with Ho = 90 km/s/Mpc on the left and decreasing 
in steps of 10 km/s/Mpc. The elliptical contour gives the 2a confidence region for the cen- 
ter of the lensing potential, as derived from radio observations of the lensed radio source by 
|Wucknitz,'Biggs, fc Browne (200 4) 

potential (see also Wucknitz' contribution to the present proceedings). The time delay 
has been measured with better than 5% accuracy (Biggs et al. 1999). 

But even this system is not 24 carat gold. The pairs of images are separated by only 
0'.'33, implying a low luminosity lens, almost certainly not an elliptical. York et al. (2004) 
have CO- added data from a large number of HST orbits (see also JACKSON AND York's 
contribution to the present proceedings). After careful PSF subtraction of the two quasar 
images, their data clearly shows an MlOl-like spiral. The residuals from PSF subtraction 
crowd the nucleus of the galaxy, making it difficult to measure its centroid accurately. 
The position of the center is crucial because of the strong dependence of the differential 
delay on the image distances. The agreement between the York et al. optical position 
and the Wucknitz center of symmetry is excellent (figure EJ. The radio position gives 
Hq = 78 ± 6 km/s/Mpc {2a uncertainty). 

The work on JVAS0218-f 357 is remarkable for the range of astronomical techniques 
and resources that have been brought to bear on the system. It was discovered (Patnaik 



Ho from time delays 
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et al. 1992) and monitored (Biggs et al. 1999) at radio wavelengths using the VLA and 
Merhn. The global VLBI data (Biggs et al. 2003) are crucial to the modeling. The redshift 

for the lens (Browne, Patnaik, Walsh, & Wilkinson 1993) required large ground based 
optical telescopes. Measuring the position of the lensing galaxy required HST (York et 
al. 2004), as is the case for most lenses. 

6. Summary and Prospects 

By the standards set forth in this paper, only one gravitationally Icnscd quasar qualifies 
as a golden lens. Taken by itself it gives a Hubble constant with a formal uncertainty 
as small or smaller than that obtained with Cepheids. But prudence demands at least 
one more golden lens before declaring that lenses do as well as (let alone better than) 
Cepheids. 

The remaining systems with time delay measurements are not without value. The cen- 
tral concentration degeneracy may reasonably be lifted by appeal to other (non-variable) 
lensed systems for which the central concentration can be measured, or alternatively, 
by appeal to measurements of gravitational potentials of nearby elliptical galaxies. The 
present author would guess that the systematic errors associated with cither of these 
assumptions would introduce an error of perhaps 5% in the Hubble constant. This would 
rehabilitate those optical quadruple systems that suffer only from the central concentra- 
tion degeneracy. 

Optical double systems suffer from both the central concentration degeneracy and the 
paucity of constraints. As they have no internal redundancy, their usefulness will depend 
upon the scatter in Hq values observed from a good-sized sample of such systems. 

In the case of B1608-I-656, the outstanding question is whether the data are so good 
that, merging lenses notwithstanding, all plausible models give the same value of Hq. 
What is needed now is the effort of a "loyal opposition" to search the far corners of 
model space. 

New telescopes and upgrades to existing telescopes are certain to produce new grav- 
itational lenses, and with them we will eventually see a Hubble constant that is less 
uncertain than that obtained from Cepheids. Even those who take it as a matter of faith 
that the universe is perfectly flat may find such a direct measurement of Hq interesting. 
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